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Beliaev damping in a superfluid is the decay of a collective excitation into two lower frequency 
collective excitations; it represents the only decay mode for a bosonic collective excitation in a 
superfluid at T = 0. The standard treatment for this decay assumes a linear spectrum, which 
in turn implies that the final state momenta must be collinear to the initial state. We extend 
this treatment, showing that the inclusion of a gradient term in the Hamiltonian yields a realistic 
spectrum for the bosonic excitations; we then derive a formula for the decay rate of such excitations, 
and show that even moderate nonlinearities in the spectrum can yield substantial deviations from the 
standard result. We apply our result to an attractive Fermi gas in the BCS-BEC crossover: here the 
low-energy bosonic collective excitations are density oscillations driven by the phase of the pairing 
order held. These collective excitations, which are gapless modes as a consequence of the Goldstone 
mechanism, have a spectrum which is well established both theoretically and experimentally, and 
whose linewidth, we show, is determined at low temperatures by the Beliaev decay mechanism. 

PACS numbers: 05.30.Fk, 67.85.Lm, 67.85.De 


I. INTRODUCTION 

The discovery of superfluidity below the A-point in liq¬ 
uid 4 He m provided a stunning demonstration of quan¬ 
tum properties of matter at a macroscopic level, paving 
the way for the experimental realization, in more recent 
times, of condensates of atomic Bose Q and Fermi [|| 
gases and of quasiparticles, like in the recently observed 
exciton-polariton condensate [|] in a semiconductor mi¬ 
crocavity. 

Even if very different in nature and in their physical 
features, these phenomena are essentially a manifestation 
of the Bose-Einstein condensation of interacting bosons 
(respectively single atoms for 4 He, Cooper pairs/bosonic 
dimers for ultracold Fermi gases, polaritons in an exciton- 
polariton condensate); thus from a theoretical point of 
view it would seem compelling to describe them, at least 
to some extent, within a common framework. 

Indeed, excitations in a superfluid can be described us¬ 
ing the quantum hydrodynamics approach developed by 
Landau [6:]; a clear advantage of this formalism is the 
possibility of describing superfluids with non-contact in¬ 
teractions and with a varying number of particles by in¬ 
troducing higher order terms by means of a perturbative 
expansion around the mean field solution. 

Collective excitations in a superfluid are destroyed 
either by Landau damping, due to their interaction 
with the thermal cloud, or by Beliaev damping, due to 
their decay into two, or more, lower energy excitations. 
There is competition between these two damping inodes: 
whereas Landau damping is relevant at finite tempera¬ 
tures, with a vanishing cross section as the temperature 
goes to zero, Beliaev damping remains the only allowed 
decay mode at T = 0. 


Therefore the Beliaev decay represents a test of Lan¬ 
dau’s hydrodynamic theory. First evidences of a phonon 
decay have been observed in superfluid liquid 4 He [2,3; 
more recently the Beliaev decay has been observed in 
a trapped Bose-Einstein condensate (BEC) of rubidium 
atoms @,E3; an analogous process has been proposed in 
order to explain the absence o f eq uilibrium in one dimen¬ 
sional interacting bosons, see [ill ] and references therein. 

In the present paper we focus on Beliaev decay and de¬ 
rive an improved version of the classical result @, UU, [l3| 
based on the observation that while the original deriva¬ 
tion requires a nonlinear term in the spectrum, nonethe¬ 
less it treats the kinematics in a low-momentum approxi¬ 
mation as if the spectrum was effectively linear. We show 
that this treatment can be extended and that, in partic¬ 
ular, the inclusion of a gradient term in the Hamiltonian 
yields a Bogoliubov-like spectrum for the bosonic exci¬ 
tations M- We calculate the decay rate for the Beliaev 
damping and show that even for low momenta and small 
nonlinearities a realistic spectrum can give appreciable 
differences with respect to the linear approximation of 
the standard result. 

The result is applied to an attractive Fermi gas: as the 
attractive interaction between atoms is tuned, the gas 
at T = 0 goes with continuity from a Bardeen-Cooper- 
Schrieffer (BCS) weakly-interacting regime, to a strongly 
interacting gas of bosonic dimers. This scenario can 
be described EM3 by introducing the complex Cooper 
pairing field, which will acquire a non-zero expectation 
value below the critical temperature. As the phase of 
the order parameter is macroscopically locked below the 
critical temperature EIG3, spontaneously breaking the 
£7(1) symmetry, its fluctuations correspond to the gapless 
mode predicted by the Goldstone theorem. These collec¬ 
tive modes turn out to be fundamental in quantitatively 
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describing the dynamics of an ultracold Fermi gas [20( | ; 
after briefly analyzing the Goldstone mode, we show that 
its linewidth gets substantially enhanced due to the Be- 
liaev decay process. We also show that our improved de¬ 
scription of the decay yields substantial deviations from 
the standard approximation. 


II. BELIAEV DAMPING: AN IMPROVED 
TREATMENT 


this constraint being fulfilled by B k = Ay. 1 . The exact 
treatment of a quantum liquid in Eq. © can be ex¬ 
panded in powers of the field operators: the first to give 
a contribution is the second order, here the theory can 
be diagonalized to a theory of non-interacting bosons, i.e. 
Hhajkbl.bk, and the requirement for H^ 2 ) to be 
diagonal fully fixes A k as: 


Ak 



h 2 k 2 \~ k 
8m 2 c 2 ) 


( 6 ) 


We briefly introduce Landau’s quantum hydrodynam¬ 
ics 0E3, a semi-phenomenological description of a su¬ 
perfluid which can be, however, rigorously justified and 
derived from the microscopical theory as discussed in [2l| . 
An exact expression for the internal energy of a classical 
liquid in a sound wave is E = f d 3 x(^pv 2 + pe), where v 
is the local velocity of the fluid, and p the local density. 
Here e represents the internal energvof the fluid for unit 
mass; Landau’s original treatment 0, [l2| assumes it to 
be dependent only on the density p , and as a consequence 
the dispersion relation for the sound waves is linear. On 
the other hand by adding a gradient term as 

h 2 (Vo) 2 

e(p)-te(p,Vp) = e(p) + A—(1) 

higher order terms appear in the dispersion relation as 
shown in jl4j], m being the mass of a fluid particle, A be¬ 
ing a dimensionless coefficient which can be fixed a pos¬ 
teriori. Within the quantum hydrodynamics framework 
the velocity and density fields of a fluid are promoted to 
quantum operators, so that the Hamiltonian for a quan¬ 
tum fluid is: 


H = 


d 3 x 


-v • pv + pe(p, Vp) 


( 2 ) 


where the term involving the velocity operator has been 
opportunely symmetrized to be Hermitean. We rewrite 
the velocity in terms of a velocity potential v = V</> and 
the density by separating the equilibrium value p from 
its fluctuations as p = p + p'. The new operators can be 
written expanding in plane waves: 

^ = 7= E ^(Ske^' + ^e-^) (3) 

IkMO 


^=iE [hB k (> eik ' r - 6 I e ' ikr ) (4) 

V |k|#0 

with V the volume of the system; the 6 q (&1) operators 
annihilate (create) a bosonic excitation over the funda¬ 
mental state of the liquid |f2), and obey the canonical 
commutation relationships. 

We impose that p' and <j> should be canonically conju¬ 
gate variables 


( 5 ) 


and the dispersion for the bosons has the usual Bogoli- 
ubov structure 


ujk = uhk 



n 2 k 2 

4 m 2 u 2 


( 7 ) 


u being the sound velocity of the sound waves in the 
quantum liquid. Clearly the original linear theory can 
be recovered by setting A = 0 and removing the gradient 
terms. The present formalism, as opposed to the Gross- 
Pitaevskii equation [H, [23|, allows for the decay of a 
collective excitation in a superfluid, in particular extend¬ 
ing the treatment to the third order one immediately sees 
that the decay of one excitation into two is allowed: this 
process is the Beliaev decay 0 described above. The 
third order term of the Hamiltonian is: 


H<A) = / d 3 r 




-A 


^- ~'\2 ^ 


8 TO 2 


(Vp') 2 ^ 


( 8 ) 


Before going on with the treatment of the Beliaev decay 
we briefly comment on the scope of application of the 
present theory; as already mentioned it can be shown 
[21| that the hydrodynamic Hamiltonian in Eq. © can 
be rigorously derived from a description of the Bose gas; 
this procedure involves integrating out the “fast fields”, 
effectively defining a momentum scale k c below which the 
perturbative expansion should be valid. Following [2lj| 
one can estimate this quantity for a weakly interacting 
Bose gas; here k c is the momentum marking the sepa¬ 
ration between a linear spectrum and the free-particle 
quadratic spectrum, and from Eq. 0 one gets 


Hk c 


2 mu 

7 T 


( 9 ) 


this condition marking, as argued in [U, the upper limit 
for the validity of the perturbation theory. 

In order to study the Beliaev decay we calculate the 
matrix element: 


= {l \ H W\ f) (io) 

between the following initial and final states: 


I*) = %|n> 

I/) = tfJlW 


[0(r),p'(r')] = -iM(r-r') 


( 11 ) 

( 12 ) 
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The matrix element in Eq. m, when Eqs. © six creation/annihilation operators; after a lengthy but 
and © are plugged in, is essentially the expectation straightforward calculation, one obtains 
value over |ff) of a number of terms composed of 

I 


h ( 3 ) = (27Tfi) 3 

fi (2V)i 


■ <$ (3) (P - Qi - 02) • 3 - qi| 


2 d u 2 


1 + Xe 3- 

dp p 


(13) 


where 9 is the angle between p and qi, the other angles 
being fixed by the condition q 2 = p — qi enforced by the 
5 function. We also defined: 


Let us make the physical meaning of the last remark 
clearer, expanding the spectrum in Eq. 0 in powers of 
k: 


Xe 1 = T ~—^—r(l + cos(0)) + cos(0) (14) 

Ip - qi| 

In deriving Eq. (fl3l) we neglected all the terms contain¬ 
ing A; it can be checked that they give oc p 7 and oc p 9 
contributions to the decay width, whereas the leading 
contribution will turn out to be oc p 5 . However the non¬ 
linear dispersion relation is relevant when discussing the 
kinematics: the differential decay rate is calculated using 
Fermi’s golden rule l34 j: 

dM , = =^\H^\ 2 S(E f Ei )-^d 3 qi d 3 q 2 (15) 
n {Ann) 


and Ef — Ei = w p — w qi — w q2 , Wk is the spectrum as de¬ 
rived in Eq. ©. The integration over d 3 q 2 is performed 
using the momentum conservation constraint appearing 
in | Hjf |. the integration over the angular part of d 3 gi 
removes the S function related to energy conservation, 
fixing at the same time the decay angle 9q, and finally 
the radial integration remains explicit. The final result 
for the decay rate is: 


w 


ii 

327 Tph 4 



qlo 


(l + Xeo^^y) 
I /'(cos 6 0 ,p,q)\ 


dq (16) 


Wk = uk + ak 3 + 0(fc 4 ) (17) 

where a has the same sign as A. The energy conservation 
constraint in the low momentum limit reads 1 — cos 9 = 
3 a{p — q) 2 and can be fulfilled only if a > 0, i.e. only 
if the spectrum grows linearly or more than linearly; for 
a < 0 no decay is allowed. 

We now focus on the strictly linear case a = 0. Energy 
and momentum can be conserved only if 9q = 0, i.e. 
the momenta of the decaying excitation and those of the 
decay products are parallel. Even for very small values 
of a the decay kinematics deviates significantly from the 
aforementioned linear situation as 9q increases. 

We stress that, even if the standard treatment of Beli- 
aev decay im correctly identifies a > 0 as a necessary 
condition for the decay to happen, then its use of a = 0 in 
the kinematics is a critical assumption; on the other hand 
the present treatment by including the gradient term in 
Eq. © allows for a realistic, Bogoliubov-like spectrum. 

Let us derive the standard result from the more gen¬ 
eral Eq. m-- having set A = 0 for a linear spectrum 
Wk = u|k| one has that 9q = 0, /' = 1 and also x = |- 
Moreover noting that / p q 2 \p— q\ 2 dq = p 5 /30, we recover 
Beliaev’s original approximation [6|, [l2| , which we report 
here for the sake of completeness: 


where |p — q|§ = |p 2 + q 2 — 2pq cos 9 0 \ for shortness sake, 
f{cos9,p,q) = (w p - w, - W| p _ q |) is essentially 

the energy conservation constraint, f is its derivative 
with respect to the first argument and 9q = 9o(p,q) is 
the only zero of / in the interval [—7r,7r], and represents 
the allowed decay angle given the incoming and outgoing 
momenta. 

Equation (1161) is the main original result of the present 
paper, which we will apply to an attractive Fermi gas. We 
stress that w in Eq. (1161) is a function of just p, u and 
of the incoming momentum p\ moreover the exact form 
of the spectrum, including the A coefficient, contributes 
indirectly to the final result, by modifying / and, conse¬ 
quently, 9q. We also note the kinematics constraints can 
be satisfied and the decay is allowed only if the afore¬ 
mentioned zero of / exists, an equivalent condition being 
that the spectrum should grow faster than linearly. 


w = p 


3207rpli 4 


1 + 


p 2 d u 2 
3 u 2 dp p 


(18) 


To conclude the present section we note that in the 
case of a weakly-interacting Bose gas Eq. m further 
simplifies, because in this case 


p 2 d u 2 
u 2 dp p 


= 0 


(19) 


Alternatively the weakly-interacting Bose gas regime can 
also be investigated, as done in jly|, starting from the 
atomic Hamiltonian, introducing the Bogoliubov approx¬ 
imation and isolating the relevant decay vertices. The 
present hydrodynamic approach is different because it 
can be derived, as already mentioned, by separating the 
fast and slow components of the fields, introducing a mo¬ 
mentum scale k c , whereas the Bogoliubov approximation 
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merely separates the zero-momentum contribution. We 
expect the two approaches to yield the same results for 
k < k c , as we verified. The hydrodynamic approach, 
however, is better suited for analyzing the collective ex¬ 
citations of an attractive Fermi gas. 


III. BELIAEV DAMPING FOR AN 
ATTRACTIVE FERMI GAS 


Let us consider a three-dimensional, uniform dilute gas 
of interacting Fermi atoms; the atoms are neutral and 
have two spin species. This system can be described 
within the path integral formalism fl5l4l7j in which the 
fermions are represented by the complex Grassman fields 
^/v(rwith the spin index er =t,4-- The par¬ 
tition function for the system at temperature T, with 
chemical potential g, is: 

Z = J ' D bl’<r>$o] eX P 5 1 > ( 20 ) 

with the following action and (Euclidean) Lagrangian 
density: 

fHP r 

S= dr / d 3 r Jzf (21) 

Jo Jv 


^ =s 22 ^° 


H 2 9 

hd T -V 2 ^ /i 

2 m 


'0a+5'V’tV’4V’4.V’t (22) 


as usual /3 = 1 /(fc^T), ks is the Boltzmann constant, V 
is the volume of the system and g < 0 is the strength of 
the contact interaction between atoms; this quartic inter¬ 
action can be decoupled, as usual, through a Hubbard- 
Stratonovich transformation in the Cooper channel, in¬ 
troducing the pairing field A(r, r) ~ Now the 

Euclidean Lagrangian density reads: 


% = 22^° 


n 2 

hd T -V 2 — fj, 

2 m 


4>a 


+A ip± + A Vr Vt - 


(23) 


In order to obtain the partition function Z one also has to 
extend the functional integration to the newly introduced 
pairing fields A, A. We rewrite the fields A, A as the 
sum of their saddle point value plus the fluctuations 


A(r, t) = A 0 + (24) 


Up to this point the theory is exact; the mean field theory 
of a Fermi gas is simply found by neglecting the fluctua¬ 
tions r/(r,r), fj(r,T). The functional integration defining 
the partition function Z can then be performed, yielding: 

Z mf = exp (-/3 fl m f) (25) 


with: 

Clmf =-^2(E sp ( k) - e k + m) - V—• (26) 

k 

with ek = h 2 k 2 / (2m)-, E sp ( k) = \J (ek — g) 2 + Aq is the 
spectrum of elementary single-particle fermionic excita¬ 
tions. The number and the gap equations for the system 
can be readily obtained from the mean-field grand po¬ 
tential 


N = - 

dg 

(27) 

dflmf' 
dA 0 , 

) = ° 

/ M 

(28) 


Lastly the gap equation in Eq. (l28l) needs to be reg¬ 
ularized, and this can be done by replacing g with the 
scattering len gth a s , according to the following prescrip¬ 
tion (see e.g. j24jb 



where a s is the s-wave scattering length associated to the 
interatomic potential. 

Let us now analyze the fluctuations contribution to the 
present theory; going back to Eq. HMD and reinstating 
the fluctuations fields r?(r, r), fj( r, r) up to the quadratic 
(Gaussian) order [23, [26| the partition function reads: 

Z = Z m f JV M ex P p«}, (30) 
having defined the Gaussian action: 

S g (v,v) = M ( ? ) ( J(- q ) ) ( 31 ) 

with q = (q, iui n ), and oj n = ^-n are the Bose Matsubara 
frequencies. The matrix in Eq. (1311) is the inverse prop¬ 
agator for the pair fluctuations of an interacting Fermi 
gas, its matrix elements are defined by [26l |27|: 


Mn(g) = ~ +W 


2 /2 
U U 


»(«) = E 


i io n — E — E' 
1 

icj n + E + E' 


2 /2 
V V 


iojn H - E E' 



1 

i uj n — E — E' 



where u = u k = ^j\( 1 + J^y), v = u k = -y/ 1 - 

u r^k+q: tj = Uk-j-q; E = E sp ( k), E = E sp ( k -{- q). 
The remaining matrix elements are defined by the rela¬ 
tions: M 22 (g) = Mn(— q), M 2 i(g) = M i2 (g). By inte¬ 
grating out the T)( r, r), fj(r, r) fields we get the Gaussian 
contribution to the grand potential [261 . 1271 ]: 

= 2)3 E ln ( det M ( ? )) ( 34 ) 
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FIG. 1: The collective mode spectrum for y = 0.0, y = 0.5, y = 1.0, from left to right. The black line represent the real part 
of the spectrum, while the dashed lines represent ± the imaginary part, the bold red line represents the threshold energy E t h■ 



FIG. 2: The pair fluctuation spectral function A vv ( k, ca) for 
y = 1, the dashed line shows the corresponding spectrum. For 
Yj > 1 the line broadening due to the Bcliaev decay effectively 
destroys a collective excitation, this is also approximately the 
scale marking the end of the validity of the perturbative ap¬ 
proach. For comparison here Eth > 3 £f 


A completely equivalent description can be given, af¬ 
ter a unitary transformation, in terms of the (linearized) 
phase and amplitude of the fluctuation field, which can 
be decomposed as ry(r, x) = (A(r, x) + i0(r, x)) /\/2. The 
Gaussian level action now reads: 


l Ed* «•) 



iM° \ 
Mft - M 12 J 


in terms of the even/odd components in i ui n 



(35) 
of the 


M matrix elements [28], [29J, i.e. M^ / 0 (q, iw n ) = 
^(M a ;,(q, iuj n ) ± M a b(q, —iw„)). This representation 
makes clear that as soon as the Cooper pairing field 
A(r,r) acquires a non-zero expectation value, i.e. un¬ 
der T c , as a consequence of the 17(1) symmetry breaking 
one expects to observe the gapless Goldstone mode [T^] . 
More specifically it can be verified from Eq. (l35l) that for 
T = 0, in the low momentum limit the phase and ampli¬ 
tude fluctuations are decoupled [29[ : the diagonal entries 
in Eq. (l35l) go to zero, and the phase (Goldstone) mode 
is gapless, while the amplitude (Higgs) mode exhibits a 
mass gap. From now on up to the end of the present 
section we will study the system at T = 0. Focusing on 
the former mode, we observe that, indeed, by solving for 
lo the equation 


det M(q, iw„ —>• w) = 0 (36) 

we obtain the spectrum of the bosonic collective mode, 
showing a gapless branch. Notably in the BEC regime 
y > 1, and across the whole crossover for low enough 
momenta, this mode takes (within very good approxima¬ 
tion) the familiar Bogoliubov-likc form 

E co i{ q) = \J e q (Ae q + 2mCj) (37) 

and the sound speed c s , along with the parameter A, de¬ 
pends on y = 1 /(fcpa 5 ). We use this spectrum in the deep 
BEC limit, while in the intermediate regime near unitar- 
ity we solve numerically Eq. (1561) to get the “exact” spec¬ 
trum within the present Gaussian approximation scheme. 
When comparing the “exact” spectrum so obtained with 
the Bogoliubov approximate form, one also has to re¬ 
member that a natural momentum scale can be defined 
by studying whether and when the dispersion enters the 
two-particle continuum reaching the threshold energy: 

-EtJq) = min(E k + £ k+q ) 
k 


(38) 
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above which a Cooper pair breaks down in two fermions. 
As far as the present work is concerned it is important 
noting that E co i grows more (less) than linearly if A > 0 
(A < 0), moreover the parameter A can be calculated eas¬ 
ily either from a numerical solution of Eq. (1361) or using 
the techniques in [30]. It turns out that A is a mono- 
tonically increasing function of y = 1 /{kpO's), where fcp 
is the Fermi momentum. In particular A takes negative 


values in the deep BCS regime and changes its sign for 
U = Vc ~ —0.14; referring to the previous section we can 
then conclude that no Beliaev decay will happen in the 
deep BCS region for y < y c . 

We now want to adapt Eq. (fl6l) to the present theory. 
We start by noticing that if the spectrum has the form in 
Eq. (1371) . then the decay angle 9 0 defined in the previous 
section has an analytic expression: 


cos 0 O (p,q) = 


m 2 c 2 

XpqH 2 


o_ 

2 P 


P_ 

2 q 


\Jm 2 cj + 2mc 2 s X (e p + e q ) - 2A pq^^/j2mc 2 + Xe p ) (2 me 2 + Xe q ) + A 2 (e 2 + e 2 ) 




139] 


We note that for the special case A = 1 , 2 mc 2 s = 2 Eq. 
(|39|) coincides with the result in 0- 
Finally the more complicated expression inside the 
parenthesis in Eq. (U6[) can be expressed using the tech¬ 
niques devised in [31| as: 

_ 30 e(y) - 8 ye'jy) - 3 y 2 e"(y) + y 3 e'"{y) 

u 2 dp p 30 e(y) - 18 ye'(y) + 3 y 2 e"{y) 

( 4 °) 

as a function of e(y) = |e F 1 £, where £ is the bulk en¬ 
ergy per particle per particle of an interacting Fermi gas; 
when calculating our final results we compared e(y) as 
fitted in m from experimental data with its mean field 
counterpart, observing no appreciable differences as far 
as the quantity Eq. (l40l) is concerned. Consistently with 
the result found in Eq. m for the weakly-interacting 
Bose gas, the quantity in Eq. 03 tends to zero in the 
deep BEC limit. 

We calculate the Beliaev decay width for the Goldstone 
collective mode of an attractive Fermi gas; as previously 
noted there is no decay in the BCS regime up to y = 
y c fts —0.14, as the spectrum as a function of |q| grows 
less than linearly. For higher values of y we can associate 
an imaginary part to the Goldstone mode spectrum, as 


hw 

Im w p = —— 


(41) 


using w from Eq. (TIBI) . In Fig. © we report the real 
spectra u> p , obtained from Eq. (15B1) . for three values of 
y = 1 /(kFds), from unitarity to the BEC regime (y = 0, 
y = 0.5, y = 1), along with their imaginary part due to 
the Beliaev decay. 

A collective excitation in a superfluid Fermi gas can¬ 
not have arbitrarily high energy, as it will be damped 
either by the dissociation mechanism at the threshold en¬ 
ergy E t h , decaying into two fermions, or by the Beliaev 
mechanism, decaying into two lower frequency collective 
excitations. Either way a natural energy cutoff can be 
associated to a bosonic excitation. 

Referring to the left pane of Fig. © we start at uni¬ 
tarity (y = 0) where the Beliaev decay width is quite 


narrow: here a collective excitation will mainly decay 
by hitting the threshold energy E t h and breaking down 
into two fermions [32 . On the other hand, approaching 
the BEC regime (y = 0.5, y = 1.0) the Beliaev decay 
width gets larger before the collective spectrum touches 
Eth- here the preferred decay mode for a collective exci¬ 
tation will be decaying into two lower frequency collective 
excitations. This trend, i.e. the progressively bigger im¬ 
portance of the Beliaev mechanism approaching the BEC 
regime, can be observed by comparing the three panes in 
Fig. ©. 

In order to define an energy cutoff due to the Beliaev 
mechanism, we can match the real and the imaginary 
part of w p similarly identifying a scale beyond which a 
collective excitation is no longer well-defined due to the 



FIG. 3: The Beliaev decay width calculated from Eq. © 
divided by the original Beliaev result (linear approximation), 
for different values of y = 1 /(fc^Os). 
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Beliaev decay. This remark is made clear by looking at 
the pair fluctuation spectral function 


A w (k,w) =-2ImG I)I) (k,w+ 7 k) (42) 

plotted in Fig. 0. As noted in j33|, it can be interpreted 
as the contribution to the density from the fluctuations 
at a given wave number q and a given momentum ui. 
In the previous equation ui is assumed to be real and 
7 k = — 7 ^ is the imaginary component of the spectrum 
due to the Beliaev decay, G vr] is the Green’s function ob¬ 
tained by inverting the matrix in Eq. m and taking 
the (1,1) entry. We observe that for low momenta most 
of the spectral weight is peaked around the dispersion 
relation, which is marked by a dashed line, assuming the 
usual Lorentzian structure. However, as the spectrum 
continues after p ~ kp, for high momenta the line broad¬ 
ening effect due to the Beliaev decay effectively destroys 
the collective excitation, and the spectral weight is dis¬ 
tributed over a large region. The border between these 
two regimes can also be approximately found by impos¬ 
ing the aforementioned condition Rewk = Imuk, which 
can be easily read from Fig. 0: when the real part of 
the dispersion is bigger than the imaginary part, the ex¬ 
pression in Eq. (1421) has a narrow peak; as the imaginary 
part of the spectrum gets bigger the Lorentzian struc¬ 
ture of the peak is lost and the excitation is no longer 
well defined. 

We can conclude that, as we go from the BCS to the 
BEC regime, the dissociation mechanism at Eth gets less 
and less relevant, as the collective mode spectrum gets 
further away from Eth', at the same time, the Beliaev 
decay channel opens at y = y c and gets progressively 
more relevant. Finally in Fig. 0 we compare the de¬ 
cay width, as predicted by the present theory, with the 
original linear approximation p, [l 2 | : even for relatively 
small momenta our treatment shows relevant differences 
with respect to the standard treatment. The differences 
get larger in the BEC regime, consistently with the fact 
that the nonlinearity term A in the spectrum is bigger; 
however we stress that even for nearly linear spectra, see 
the cases y = 0 and y = 0.5 in Fig. 0, the correction 
due to the present treatment can amount up to 25% for 

ifcF - 1 ' 

In conclusion we briefly comment on the scope of ap¬ 
plicability of the present theory to the fermionic case; 
adapting Eq. 0 one finds 


k c — 


2 mc s 


(43) 


and we do not expect the theory to be applicable above 
this momentum threshold; a direct calculation shows that 
starting at unitarity, up to the moderate BEC regime we 
considered in Fig. 0 and Fig. 0 , k c assumes respec¬ 
tively the following values: k c = 3.0 6k F , k c = 1.72 kp, 
k c = 1.31fc F ; in the deep BEC limit we get 

fc ~ (44) 

c ~ 2 m (44) 


We notice that for the cases we considered, the momen¬ 
tum scale k c marking the breakdown of the perturbation 
theory is higher or equal to the scale at which a collec¬ 
tive excitation is no longer well defined due to a high 
decay rate; we conclude then that the present treatment 
is consistent. 


IV. CONCLUSIONS 


We have extended the standard result of Landau’s 
hydrodynamic theory of a superfluid, which leads to a 
purely linear spectrum implying a collinear Beliaev de¬ 
cay. By including a gradient term in the Hamiltonian, 
we have recovered the Bogoliubov-like spectrum bosonic 
excitations in a superfluid have a Bogoliubov-like spec¬ 
trum of excitations producing a larger phase-space for 
the Beliaev. We have shown that even slight variations 
from linearity of the spectrum can give important modi¬ 
fications to the decay rate of the process we consider. 

We have applied our result to an interacting Fermi gas 
in the BCS-BEC crossover: we have shown that no de¬ 
cay happens at zero temperature in the deep BCS regime, 
due to kinematics constraints, as the spectrum grows less 
than linearly. As the strength of the attractive interac¬ 
tion is increased, the collective mode spectrum increases 
linearly or faster as y > —0.14, thus allowing the decay 
of one collective excitation into lower energy excitations 
and this mechanism becomes more and more relevant as 
the coupling gets stronger. We observe that in the BCS 
regime in the low-temperature limit a collective excita¬ 
tion can decay only by breaking down into two fermions 
at E t h', on the other hand at unitary and in the BEC 
regime a collective excitation can also decay in two col¬ 
lective excitations by means of the Beliaev decay. 

Finally we have identified the regimes to which the the¬ 
ory we have developed applies. The perturbation theory 
behind the hydrodynamic treatment of a quantum liq¬ 
uid breaks down at a critical momentum k c . We have 
estimated this value, verifiying the internal consistence 
of our treatment, showing by analyzing the decay width 
and the pair fluctuation spectral function that in the 
fermionic case k c is higher or equal to the momentum 
scale at which a collective excitation is no longer well 
defined due to the decay process. 
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